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ABSTRACT 

The diffusion equation for cosmic ray nucleons in the 
atmosphere has been solved analytically, taking into 
account the transverse momentum distribution of nucleons 
produced in nucleon-air nucleus collisions. The effect of 
the transverse momentum distribution increases the nucleon 
intensity at large zenith angles and low energies. 


1 . Introduction 

The intensity of cosmic ray nucleons in the atmosphere, due 
to cosmic ray primary protons, can be expressed by (1) 

N ( E , x ) = C p E _(Y +1) e _x/A (1) 

where N(E,x = 0) = C E + is the primary proton spectrum, 
x the atmospheric depth and X the diffusion length for 
nucleons in the atmosphere. This solution is obtained by 
the assumption that the secondary nucleons are moving in the 
same direction as their parent nucleons.lt can be used for 
deriving the intensities of cosmic ray mesons in the 
atmosphere and muons at sea— level ( 1 ) because the energies 
of interest here are high such that the average transverse 
momentum of secondary nucleons <p > is very small compared 
to the energy E of the nucleon, <pS << e. If <p > is no 
longer much smaller than E, as is the case, for Instance, for 
16w energy cosmic rays in the upper part of the atmosphere, 
equation (1) will give too low nucleon intensity at large 
zenith angles 0, as we shall see below. 


(3.) we obtain the 


2 . The Atmospheric Depth 
From the U.S. Standard Atmosphere (197£) 
atmospheric depth at sea-level and zenith angle 0°, 
x o = 1035,6 9 /cm . We approximate the Standard Atmosphere by 
an exponential atmosphere and obtain for the atmospheric 
depth at hight Z (km) above sea-level x = x Q exp(-z/H), 

H = 6.95 km. The error is less than 12% for x < 1 g/cm* 
compared to the Standard Atmosphere. For an inclined 
direction, 0° < 0 < 90°, we obtain for the atmospheric depth 


x m ax ( 0 ) / at sea-level. 
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p Q is the density of air at sea-level 
R is the radius of the earth, 6366 km. 

The error, due to approximations, is less than 5%. 
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From the figure we obtain: 

cos ( 0 ) = cos ( 0 ) cos(dJ-) + sin (0) sin^?-) cos(<^>) 
will be given by later c 

hi 

and is small so that we can take 

a <Lp t>© ->■ 9 o>~ Co y 

We introduce the reduced zenith 
angle S, 

« "if/ - 

1 \1 + 4 tfe 06 *] 


which is the local zenith angle at 

approximately half the total atmospheric 

depth x (0) . 
max 

Let x be the atmospheric depth at hight Z 
above sea-level amJ zenith angle 0. Then the 
atmospheric depth x at hight Z and zenith 
angle % is given by f - f 
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3. The Diffusion Equation 

Let the differential production cross section for nucleon 
production by nucleons have the following transverse momentum 
distribution £ — <* z. 

F't — =- (L ’ cA ■= CcTWo'fc . 

where p„ is the longitudinal momentum and p is the trans- 
verse momentum of the produced nucleons. With A>-= we then 
obtain for cosmic ray nucleons g; / E \2. 

with g nn = j <$C6 ,£'<*) <2^ t 't; 

O 

We use the generation method and obtain the diffusion equat- 
ion for the first generation nucleon 


4f--x; IS e "V £ ' e> ' t> ■ 
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. H Z J = 84 g/cm 2 is the inter- 
ior f c\ ( lt- g, action length for nucleons. 
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The analytic solution is 
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The solution for the m'th generation nucleon has the form 
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If we define an average for g , for instance: 

3 m,n 

we can sum over m and write an approximate solution for the 
nucleon intensity in the atmosphere 

A/ Cb,x,6) ^ c e e ’ 


and i^ t y e O) i s i ncreas i n 9 with 

' S l/*l I V 

-ck^L. 


increasing and zenith angle &. 

-V £= 

If the differential production cross section is rV €■ 
we obtain for the first generation nucleon 
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f ; % y. 

and we see that this flatter trans verse momentum 
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distribution gives a greater effect. 

4. Discussion and conclusion 
" — — 

N(E,x,0) is the solution of the 3 - dimentional diffusion 
equation for nucieons in the atmosphere and N (E,x) the 1- 
dimentional solution^ as given above. The figure shows, as an 
example, the ratio N (E,x,&) /N (E,x) at 30 km above sea-level 
(x ( & = 0° ) = 12.4 g/cm^) for & = 90° , x = 350 g/cm^ as a 
function of energy E, 



We see that below 10 GeV the ratio ^/N ij^ increasing rapidely 

with decreasing energy. This is because N/N is a function of 

“ 2 ^ 

E as we see from the solution for N' above. 

In the case of ^ , we obtain approximately 

the same curve, but with the energy a factor 2 higher and 
the rapide increase is between 20 GeV and 10 GeV. We see that 
the effect of the transverse momentum distribution is not 
only a function of the average transverse momentum <p >, but 
depends also on the shape of the transverse momentum 
distribution . 
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